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1. Introduction 

Fermions coupled to inhomogeneous background bosonic fields can become localised, in 
the sense that the fermion spectrum contains localised bound states [Q, ^. In particular 
in the presence of scalar kinks, fermions prefer "living on" the kinks rather than populate 
higher energy (delocalised) radiation modes. This phenomenon is ubiquitous in condensed 
matter systems with impurities and where external magnetic (gauge) fields provide the 
inhomogeneities. The case of scalar fields acquires particular interest, if one allows kinks 
to represent higher dimensional domain walls or branes in string theory [Q]. The localisation 
of fermions on inhomogeneities then becomes reminiscent of brane world scenarios, where 
the Standard Model fields are expected to "live" on one particular three dimensional brane. 

In the presence of extra dimensions, additional branes may exist (in the "bulk"), each 
with their localised fermion states. As such branes collide, one may envisage fermions being 
transfered from one brane to another. Such a scenario was studied in |p. 

Another picture of braneworlds emerges from the work of Horava and Witten , where 
it was realized that spacetime boundaries can play an important role. In that particular 
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case it was discovered that each of the two boundaries supported an E8 gauge theory. This 
braneworld model then places our Universe at the boundary of some larger spacetime. The 
question of what happens as bulk branes collide with our boundary-Universe is what this 
paper addresses, as well as extending the calculation of Gibbons 

In order to model a 3+1 Universe as a boundary of a 4+1 spacetime we could start 
with an action of Dirac fermions, and scalar fields that support domain walls to model the 
branes. However, if we employ a planar symmetry along the walls this model reduces to a 
1+1 model with a boundary. And for our particular initial conditions (domain walls with 
just a bound zero mode) the equations for the Dirac fermion reduce to the equations for 
a Majorana fermion. So, while we in practise simulate a Majorana fermion in 1+1, this 
can be lifted to a 4+1 spacetime with Dirac fermions, a scalar field for the branes, and a 
boundary. 

We study the fermion transfer numerically, treating the scalar as classical and the 
fermions in terms of a set of quantum modes. 



2. Scalar-fermion model in 1+1 dimensions 

We consider a model of a real scalar (j) and a single fermion species ip in 1+1 dimensional 
space-time. In our simulations we shall be considering situations both with and without a 
boundary, in either case the bulk action we use is 

^1„ ,.-„„, A . 

4^^ 



'S'buik = - J dtdz 
with the bosonic boundary action being 

'S'boundary — -f I dt 

Our choice of conventions is 



1 



(2.1) 



2 = 



(2.2) 



(2.3) 



At this point, the couplings A and g are free and we have chosen the boundary to lie 
at z = 0. The boundary action for the scalar may seem unnatural at first but, as we 
shall see, such a term means that there is no force between the boundary and a (static) 
kink, (or antikink if the + sign is chosen). Another reason is that this bosonic boundary 
term is required for the action to be supersymmetric, along with taking the spinor to be 
Majorana with couplings related by g^ = 2A. We shall not be requiring supersymmetry. 
The corresponding equations of motion in the bulk are 



i7''d^ + gcf>iz,t))i^iz,t) = 0, 
[d^d>' - A {cf>\z,t) - 1)] (t>{z,t) = -igi,^{z,t), 



By choosing a real representation of 7^, 



70 





(2.4) 
(2.5) 



(2.6) 
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the equation of motion for the complex two-component fermion sphts up into two uncoupled 
real (Majorana) copies, 

V^ = < + iV'f. (2.7) 

The equations of motion for il^i^ and are identical, and we will from now on think in 
terms of a single, two-component Majorana fermion -0^^, ignore the superscript, and write 
it as 

<(z,t)=f^^5"'';+^^5"'';V (2.8) 



The equations of motion then read ^ 



Mz,t) = -dMz,t) + gU2{z,t), (2.9) 
Mz,t) = -d,Mz,t)- gHi{z,t), (2.10) 
0(z, t) = dl<\>{z, t)-X (02(z, t) - 1) ^{z, t) (2.11) 

We discretise these in a straightforward way and solve them numerically using a standard 
algorithm, 4th. order accurate in time, 2nd order in space. The boundary condition was 
satisfied using a Newton- Raphson iteration [|l^] . 

The boundary conditions coming from the action are 

9.</'|o = Ty|(0'-l)|o, (2.12) 
and for the fermions we will be using 

tpilo = Ttpilo, tp2\o = ±V'2|o, (2.13) 

which can be derived from a boundary action of it | J dt['i/'V']z=o- We shall call these the 
=p Boundary Conditions, ^BC. We see therefore that at the boundary either ■01 or ■ip2 
vanishes, depending on the choice of sign of boundary conditions, with the —BC giving 
"0110 = and the +BC giving -0210 = 0- As we shall see later, these boundary conditions 
allow for a normalizable fermion condensate on the boundary to co-exist with a kink or 
antikink. 

The scalar field equation can be thought of as concerning a classical field, but as 
fermions are quantum, the right hand side source term should be replaced by —ig{'4)il)). 
The resulting equations amount to a quantum fermion in a classical scalar background. 

For the purpose of this paper, we will make the approximation of neglecting the fermion 
back-reaction on the scalar. The upshot of this approximation is that whereas the solution 
of the mode equations (see below) are independent of the initial state (i.e. particle content) 
of those modes, the back-reaction term (V'V') is a quantum average over some density 

^ These equations are equivalent to those of the 4+1 system in if we take ipi = iip^ and ^2 ~ ip-- 
This explains why the authors see an additional phase of 7r/2 between tf}+ and tp-. 
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matrix. In the present context, we are more interested in the behaviour of the fermion 
modes as the scalar background changes than the exact details of the kink evolution. It 
is however clear that with many fermions present, in particular when including all the 
non-localised modes, the back-reaction may be sizeable, and may even drive the system to 
a high-temperature thermal state. Then the setup of a solitary kink or kink-antikink pair 
may no longer be realistic. 

In our simulations we monitor the conservation of energy and we give here the expres- 
sion of the energy in the scalar field. 



H = dx 



±W{z = {)). (2.14) 



In particular we see that the boundaries have an energy associated with them. In the case 
where the scalar field is in the vacuum at the boundary we see that the boundary has an 
associated energy of 

Sf,(0 = l)=T^yA72, ^fe(0 = -l) =±^072- (2-15) 

This boundary energy will be important in understanding the dynamics of kink-boundary 
collisions. 

2.1 Kinks and boundaries 

The static scalar equation of motion has a kink and an antikink solution (pK and (j)A, 

<Ai,(z)=tanh(^^), 0^(z) = -tanh(^^), = ^f, (2-16) 

where zq is the center of the (anti-)kink. We will use A = 2, L> = 1 throughout.2 Note that 
these static solutions obey 

9.0 = Ty|(<^'-l), (2.17) 

so that there is no force between a kink and a —BC. Similarly there is no force between 
an antikink and the +BC. In the case where the coupling constants take the supersym- 
metric values, = 2A, the (anti)kink and the {-\-BC)—BC break the same half of the 
supersymmetry. 

When colliding kink with antikink we shall employ periodic boundary conditions. 
When colliding single kinks onto a boundary, these are placed at z = 0, with the kink 
coming in from the left {z < 0) . Note that there is no loss of generality by choosing to send 
in kinks but not antikinks, while using both boundary conditions. The system of kink and 
—BC is equivalent to antikink and +-BC, while kink and +BC is equivalent to antikink 
and —BC; this covers all the possibilities. 



This amounts to rescaling the original action (2.1) 
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3. Bound states 



3.1 Static kinks 

An interesting facet of topological defects such as kinks is that they often allow for bound 
states in the particle spectrum. In the case at hand there are bound states for both the 
scalar and fermi fields. Indeed, the existence of the scalar bound state was the motivation 



behind the work of Rubakov and Shaposhnikov |13], asking whether we live on a domain 
wall. To study the scalar spectrum we consider perturbing the scalar equation of motion 
in ( |2.5| ) around a kink 

(p{z,t) = (/)kink(2) + (5(/>(z,t), (3.1) 



then by writing 



S4) = exp{iujt)F{z), (3.2) 



we solve the resulting eigenvalue equation to find the frequencies of the bound states [p^, p^. 
Taking the change of variables ||6| 

Z = t&nh{z/D), (3.3) 
we arrive at an associated Legendre equation with 1 = 2, 

for which there are three solutions, P\{Z), P\{Z), P%{Z), corresponding to m? = 4 — 
uP'D'^ = 2, 1, 0. These give frequencies w = 0, ^/?,/D, 2/D. The first of these modes, 
the zero mode, corresponds to the translation of the kink, the second is a true bound state 
while the third is not normalizable and so is not considered to be in the physical spectrum. 
We also note that this last solution is right on the boundary of the continuum states, which 
have mass \/2A = 2/D. 

As well as these scalar bound states there are fermion modes localized to the kink 
which we now describe. First there is the zero-mode Ol 



V;f(z,0) = J ^[^-^ + y^l \, , i^^{z,Q)=Q, (3.5) 

where we have normalized too unity according to the inner product 

(V',X) = / dx'4>'^x- (3-6) 



Now recall that the scalar field of the kink can happily co-exist with the -BC, and that 
these boundary conditions required the vanishing of V'2(-z = 0). So we see that the -BC 
also have no effect on the fermion fields of the static kink. 
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In addition to the zero mode we also find that there is an excited fermion mode given 

by 



ijf'^{z,t) = -MujD- ;_^J^j^ cos{ujt-^), (3.7) 



[cosh(^)]' 



V'f ^(z, t) = N- ; sinM - (3.8) 



1 



[cosh (^)] 

with 



2 _ 2gZ?-l 2 _ {gD-l)T[gD+ 1/2] 



Recall that this model, with a Majorana fermion, is supersymmetric when g = v2A = 2/1). 
In this case we see that the frequency of the excited fermion mode matches that of the scalar 
bound state, uj = V^/D, as is expected for super symmetry. We also expect, given that 
there are no more scalar bound states, that there will be no more fermion excited states. 
We note that this does not hold true as we change couplings away from the supersymmetric 
case.^ 

3.2 Moving kinks 

As we are interested in kinks colliding with each other and against boundaries we need to 
know what kinks and bound states look like when they are moving, that is, we need to 
boost the kink. 

The scalar field of a kink (antikink) moving at speed v simply involves a Lorenz con- 
traction by 7 = l/Vl — v'^, 



(^f/^(z,0) = ±tanh 



An isolated, moving kink or antikink has an energy of 



(3.10) 



4 

-E^antikink = -E'kink = ^ V-^/^T- (3-11) 

We note that the energy of a static kink, v = coincides with the difference in energy 
between a — and a + boundary, AE\^ = | y^A/2, from Eq. ( p.l5| ). 

The fermion modes also transform in the usual way, and in terms of the 1, 2 components 
we have 

Mz,t)^i^l,iz,t) = ^|2±l(^^,[J{z-vt)] + -^^P2biz-vt)]y (3.12) 
Mz,t)^^P'^^,iz,t) = ^2±l(^iJ^[J{z-vt)] + -^^Pl[^iz-vt)]y (3.13) 



^We would like to thank Kei-ichi Maeda for informing us of his calculation showing that there is a tower 
of bound states with frequency lj^ — n ^^^^" , for integer n subject to < n < gD. 
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3.3 Boundaries 



Just as kinks support bound states, so does a boundary. To describe these we consider the 
scalar field to be in one of its vacuum states at the boundary, (j){z = 0) = ±1, then we see 
from the spinor equations of ( |2.9| ) that the boundary carries a single localized, normalisable 
fermion 

d,Mz) = -gct>{z)iJi{z), dMz) = g<P{z)Mz)- (3-14) 

For (j){z = 0) = +1 the normalisable solution is given by 

^Pfiz) = 0, ^P^{z) = V^exp {gz) . (3.15) 
For (j){z = 0) = — 1 we instead have 

V'f (z) = V^exp {gz) , ^l;^{z) = 0. (3.16) 

Note that the condition of normalisability imposes that one fermion component is zero in 
each case. For example, if we have a single kink in the bulk (which allows a condensate of 
V'l) then at the boundary we will have (/> ~ +1 which allows a boundary condensate of xp2, 
consistent with the -BC. 



4. Kink dynamics 




Before discussing the behaviour of fermions 
on kinks we first want to describe how kinks in- 
teract with antikinks, and also how they inter- 
act with the zizBC boundaries. The dynamics 
of kink-antikink collisions have been studied by 
a number of authors [0, 10, revealing 

the rich structure of behaviour at small impact 
speeds, v. At large speeds, v > 0.2, the kinks 
simply collide a single time and bounce away to 
infinity. As we go to smaller speeds then the 
kinks may have multiple collisions before trav- 
elling off to infinity, or they may simply annihi- 
late. While extending the analysis of Q for the 
fermions involved in such collisions we shall con- 
centrate on the simplest range of speeds where 
there is a single collision. Fig. ||. 

In order to capture the dynamics of brane-boundary collisions, Antunes et al |16| 
modelled the system with a scalar field and, in our language, collided kinks against zizBC- 
like boundaries. They discovered that the kink is temporarily absorbed into the boundary, 
but re-emerges having lost a (model- independent) fraction of its kinetic energy of ~ 63%. 



Figure 1: The incoming kink profiles be- 
fore, during and after a kink-antikink colli- 
sion. The curves are equally spaced in time, 
with dashed lines before the collision and 
full lines after. 
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Figure 2: The incoming kink profiles before and after tfie collision with a —BC boundary (left). 
The curves are equally spaced in time, dashed before the collision, full after. The kink enters the 
boundary and comes back out unscathed. In contrast, the +BC boundary (right) decays straight 
away and emits an antikink which prevents the incoming kink from reaching the boundary. 

The boundary conditions used were in fact slightly different from ours, in that they made 
the replacement 



hence having —BC between the two potential minima and +BC outside. This fix is 
responsible for the reported energy loss at collision. 

In our simulations we shall be re- visiting this scenario, but using strict zizBC bound- 
aries. Then energy is conserved at collision, taking into account changes in boundary 
energy. For —BC, the kink reemerging from the boundary has equal and opposite velocity 
to the incoming one. Fig. Q (left). In the case of +BC the vev (/>(z = 0) = 1 at the bound- 
ary is unstable to decay via the emission of an antikink, Fig. § (right). This because the 
energy difference between +1 and —1 on the boundary is exactly the energy of an antikink. 
In the absence of an incoming kink (or with a kink very far away), the boundary could 
remain in the unstable vacuum, but as the kink approaches, the exponential tail hits the 
wall and causes the vacuum to decay. In practice, this means that before the kink reaches 
the boundary, an antikink will be emitted and collide with the incoming kink. Hence we 
cannot realise a kink/-|-i?C collision. 

5. Particle number and Bogoliubov coefficients 
5.1 Kink-antikink modes 

The general philosophy to calculating particle numbers is to identify the correct vacuum 
and creation/annihilation operators. For example, when colliding a kink and antikink we 
consider the system in the asymptotic past to be in a vacuum state and we can therefore 
expand the fermi wave operator as 





* = a'^i^g + a^'^^g^ + a^i^t + a'-^i^f^ + continuum. 
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Where the a are the particle operators for: the fermi zero mode on the kink (K); fermi zero 
mode on the antikink (A); first excited fermi mode on the kink (KE); first excited fermi 
mode on the antikink (AE). And the ipin are the (normalized) mode functions found in 
section ^ corresponding to kinks and antikinks with the requisite position and velocity. If 
there is more than one excited fermion mode then they may also be included in an obvious 
way. To an observer in the asymptotic future there will be a similar expansion, only now 
they will use a different set of particle operators, 

^ = b''^b!^ut + b'^'^^P^J + 6^Volt + b^^i^^^t + continuum. (5.2) 

Again, the mode functions ipout are the mode functions corresponding to the kink, antikink 
with the appropriate position and velocity of the outgoing defects. 

From the original action we see that the momentum conjugate to the wave operator 
^ is so by using the standard equal-time anti-commutation relation 

{'i>a{t,x),^l{t,y)} = 6a,(i6{x- y), (5.3) 

we see that the particle operators obey 

{a,a^ = {b,b^ = l, (5.4) 

with other anti-commutators vanishing. 

Now that we have our wave operator in the asymptotic limits, we need to relate the 
particle operators in order to understand how particle numbers are affected. To do this we 
introduce the Bogoliubov coefficients in the standard way, which relate the mode functions 
ipin to ipout in the asymptotic future. In the following expressions the mode functions ijjin 
are the time-evolved mode functions from ( |5.1| ) evaluated in the asymptotic future, 

i^fn = aKi^out + Mouf + iK^out + ^Ki^out + continuum, (5.5) 

i'i^^ = OiKE^out + PKElPout + iKEi'out + ^KElPout + COntinUUm, 

i^fn = OiAi^out + ^Ai^oJt + lA^^t + ^A-^^t + continuum, 
i^tn = otAE'4'out + f^AEi^ouf + lAE'^out + ^AE^i^t + continuum, 

The Oj, 7j, 5i are the Bogoliubov coefficients, which can be extracted by taking the 
inner products according to (|3.6D 

and similarly for ipKE, i^A, i^AE, i^B- In this way we are able to calculate all of the 
Bogoliubov coefficients for a given simulation. 

The number operator takes the standard form, 

iV = i / da;(^'^^' - ^'^'t), (5.8) 
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which we may write, using ( |5.lD , (5.2), (|5.4| ), in terms of the particle operators 



n 



AT = +n^^ + n^ + n^^, (5.9) 

f^KE ^ ^Km^jKE _ 1 ^ ^AE ^ ^AE^ ^AE _ (5 -^^) 

where O represents either the a or 6 particle operator depending on whether we are looking 
at the past or future respectively. We define the vacuum relative to the initial state as 

Vi,j, a^|0) = 0, (5.12) 

and the particle states as 

li^OOO) = a^t|o), loAOO) = a^t|0), |00i^£;0) = a^^^\Q), \m^AE) = a^-^t|o). (5.13) 

We can see from this that, as described in [P, the vacuum states have fermion number — | 
whilst the excited states have fermion number 

We are now in a position to express the asymptotic-future particle operators, 6, in 
terms of the asymptotic past particle operators, a, by comparing (5.1) and ( p.2| ) and using 
(O), 



= axa^ + aKEa^^ + otAa"^ + aAEa"^^ , (5.14) 

h^^ = /3^a^ + /3x£;a^^ + /J^a^ + ^AEa"^"", (5.15) 

= jKa"" + ^KEa""" + 7Aa^ + 7Ai?a^^, (5.16) 

6^-^ = Sxa^ + 6KEa^^ + (^^a^ + 6AEa^^ ■ (5.17) 

For simplicity, we will focus our interest on evaluating the expectation values in the 
state corresponding to a kink with a zero mode fermion colliding with a vacuum antikink. 

= {KOOO\b^^b^\KOOO) - 1/2, (5.18) 

and similarly for n^^, n^, n^^, . These then represent the fermion occupation num- 
bers in the respective modes after the collisions, assuming that the initial state has only 
excitations (particles) in the kink fermion groundstate. This state is annihilated by all the 
^KE/A/AE ^ leaving us to calculate only 



n^ + 1/2 = (n^ + 1/2) |a-^|2, (5.19) 
n^^ + 1/2 = [4 + 1/2) (5.20) 



n 



^ + 1/2= (n^ + l/2)|7^|2, (5.21) 



n^^ + 1/2 = (< +1/2) (5.22) 

where we have generalised to an initial state with fermions only in the K mode, but with 
an arbitrary particle number Uq . We will not be concerned with the normalisation of the 
state, but simply compute the Bogoliubov coefficients ax, Pk, Ik-, 5k- Below we will 
suppress the label K. 
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5.2 Kink-boundary modes 

For the collision of a kink on a —BC, we now expand the wave operator as 

^ = a^i^^ + a^^t/;^^ + a^V»^ + continuum. (5.23) 

and 

^ = b^'i^^^t + b'^^'tl^^J + b^'i^^^t + continuum, (5.24) 

where the index B refers to the boundary zero mode. The expansion of the mode functions 
proceeds as before, with the addition of a boundary mode function, 

= aK'4'out + Pxi^aaf + ^Kipout + coutinuum, (5.25) 
i'Fn = o^BtPout + (^BiPout + ^Bipout + contiuuum, (5.26) 

and if we send in a kink with non-vanishing occupation number for the zero-mode 
fermion, then we find that the boundary number operator in the asymptotic future is 

n^ + 1/2 = (n^ + 1/2) (5.27) 
6. Kink-antikink collisions 




Figure 3: Left: The scalar field (black) and the fermion K mode components tpi (red) and ijj2 
(blue) before (dashed) and after (full) the collision. The initial velocity was v = 0.6, and the 
coupling 5 = 2. Right: The overlap of the incoming K mode after the collision with the outgoing 
kink K and antikink A modes (a and 7), and the KE and AK modes (/? and 5). 

The first set of results that we shall present is an extension of the study in j^, where 
the Bogoliubov coefficients of the fermion zero mode were calculated for kink/antikink col- 
lisions. Here we also present data for the excited fermion bound state, as well as observing 
the dependence on collision speed. 

The collisions were performed by initially placing a kink and a antikink a distance 30v 
apart. We then boosted them with velocity v and —v respectively, as described above. 
Fig. m (left) shows the profile of the scalar field and fermion modes before and after the 
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collision, for the case of f = 0.6, g = 2. Both before and after the collision, the fermion 
modes are well localised around the kink and antikink. 

We then calculate the Bogoliubov coefficients in time, Fig. |3| (right). Before the 
collision at Dt ~ 15 the fermion mode is the kink K mode, and so |ap = 1. during the 
collision all bets are off; in particular we are not able to assign velocities to individual kinks. 
Already at time Dt = 20, the kink-antikink pair have disentangled themselves, and a final 
Bogoliubov coefficient had been established. Although there is some residual oscillation 
even at late times, we assign final values at Dt = 50. 

It is worth noting, that because the phase f in Eqs. {^^, 3.8) is undetermined, we 



instead keep the whole time-dependence Lvt + cp fixed when computing the overlap. This 
means that the computed /3 and 5 are oscillating functions in time. Fig. |3| (right), and 
we should use the amplitude of this oscillation as the Bogoliubov coefficient. Taking this 




g 



Figure 4: Left: The velocity and coupling dependence of the Bogoliubov coefficients of the incom- 
ing K mode with the outgoing K mode. Right: The same for the overlap on the A mode. Errorbars 
reflecting the residual oscillation (see text) are roughly the size of the symbols. 



into account, the Bogoliubov coefficients can be read off with an accuracy of in most cases 
better than 0.01. Fig. Q shows the g and v dependence of and 7^, the overlap with the 
outgoing K and A modes. 

As noted in [^, the coupling dependence approximately follows a a-\-bsin{cg + d) form, 
although the amplitude decreases somewhat with g, especially for a. We will not attempt to 
fit this behaviour numerically, but simply note some qualitative points of interest, and 
7^ are anti-correlated with the same period in g, and the bulk of the fermion number ends 
up in these lowest energy modes after the collision. However, as v is increased, fermion 
number is lost from these modes, in particular from a^. For small g, the Bogoliubov 
coefficients take a very long time to settle and seem to decrease continuously until they 
do so. Fig. P shows the corresponding overlap with the KE and AE modes, and 6"^. 
These are strongly correlated with each other and seem to have the same period in g as 
K/A mode coefficients, but with a phase shift of tt/2. The fermion number taken away in 
these modes is much smaller, but it is interesting that up to 10 to 20 percent can end up 
here. There is a mild increase with v consistent with the decrease in and 7^. 
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Any missing fermion number must ther 
be transferee! to modes which we do not 
take into account, radiation or additional 
time-dependent bound states. We quan- 
tify this by calculating the sum of the co- 
efficients |ap + -I- I7P -I- shown in 
Fig. m. The loss to radiation has some non- 
trivial dependence on g, and increases with 
V, presumable because there is then energy 
available to excite higher energy modes. It 
is striking that up to 60 percent of the 
fermion number can be lost to radiation 
in this way. 

7. Kink-boundary collision 

The second set of results in this paper are those associated with a brane colliding into a 
boundary, rather than another brane. as such, we now collide a single kink onto —BC 
boundary and see to what extent a fermion originally localised on the kink will stick to the 
boundary. Fig. ^ (left) shows the field profiles initially and at late times after the collision. 
Again, the fermions remain nicely localised around the kink, but after the collision also at 
the boundary. 

We shall again consider only the case where the kink starts with a fermion zero mode, 
and calculate the Bogoliubov coefficient for finding a fermion in the outgoing K (|ap), the 
outgoing KE and the boundary B {\^\'^) modes. Fig. ^ (right) show the coefficients 

in time, again with = 1 until the collision at Dt ~ 15. After another transient collision 
stage, the coefficients take on definite values, and we end the simulation at Dt = 50. 
The results concerning fermions radiating into the bulk are somewhat different to the 
kink/antikink collisions. As shown in Fig. ^ we see that the sum representing bound state 
fermions, |ap + |/5p-|-|^p, is very close to unity. This implies that very little of the fermions 
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Figure 6: The sum of Bogoliubov overlaps on 
all the bound states jap + |/3|2 + j^p + |J|2, The 
deviation from 1 is the amount of fermion number 
carried away as radiation. 
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end up in the bulk, with larger collision speeds producing more bulk fermions as one would 
expect. We again determine the dependence on g and v, shown in Fig. |8|. In this case is 
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Figure 7: Left: The scalar field (black) and the fermion K mode components ipi (red) and -02 
(blue) before (dashed) and after (full) the collision. The initial velocity was v — 0.6, and the 
coupling g = 2. Note the non-zero boundary mode contribution to the far right. Right: The 
overlap of the K mode after the collision with the outgoing kink K and KE modes {a and f3) and 
the boundary mode B (^). 

almost exactly anti-correlated with Most fermion number is transfered to the boundary 
for small and large values of g, with the maxima moving down and up , respectively as v is 
increased. For v = 0.9, the kink retains its fermion, at least in the range of couplings used 
here. 




Figure 8: The Bogoliubov coefficients of the incoming K mode unto the outgoing K mode, |ap, 
(left) and B mode, (right). 



In order to understand the kink/antikink results we repeat the approximation devel- 
oped in Q. The approximation is a way of solving ( |2.5D , ( p. 10 ) which writes the fermi field 
during the collisions as 

^^(1,2) ^ Al,2)/(^). (7-1) 

where f{z) is some even, normalized function, and the total amplitude is normalized by 
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A^-^^ + A^^) — 1- We then find that ( |2.5| ), (|2.10 ) may be integrated along the z axis to give 

^(1,2) - ^ (1 + sin(2g0eAi)) . (7.2) 

In this expression we have a representative value for during the collision, (j)c, and a collision 
timescale At. In this way Gibbons et al were able to explain the sinusoidal behaviour of 
the Bogoliubov coefficients when a kink collides with an antikink. Unfortunately, the case 
of a kink colliding with a boundary does not succumb to the same analysis, largely because 
the fermion mode functions have rather different z-dependence owing to one of them being 
forced to vanish at the boundary. However, from the form of the Bogoliubov coefficients it 
is tempting to speculate that a relation similar to ( |7.2[ ) holds. A possible explanation for the 
longer "wavelengths" in Fig. |8| could then be that during the kink/boundary collision the 
value of (p changes very little from its vacuum value, while in the kink/antikink collisions 
one finds that <j) overshoots the vacuum by some amount depending on collision speed. 
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Figure 9: The sum, jap + |/3 
ing the fraction of fermions that end up in bound 
states. 



8. Conclusion 

In summary, we have performed a detailed 
numerical study of fermion transfer in kink- 
antikink and kink-boundary collisions. 

In kink-antikink collisions, we confirm 
the findings of Gibbons et al Q that al- 
though the scalar field kinks bounce off 
each other in an elastic way, fermions ini- 
tially in the K mode on the kink will be dis- 
tributed on the K, A, and radiation modes. 
As an extension of their work we also in- 
cluded the first fermion excited modes, KE, 
AE, finding that these modes also gets excited, and presented data for a wide range of colli- 
sion speeds. The distribution is very sensitive to the value of the coupling g, and hence the 
mass of the fermions g(j). For small incident velocities v, most of the fermion number ends 
up in the K and A modes, and these are anti-correlated in g. As the velocity is increased, 
more and more fermion number is transfered from the K/A modes to the KE and AE, but 
also delocalised radiation modes, and up to 60 percent can be "lost" to the bulk. 

When colliding kinks on a -EG boundary, we found again that the kink is reflected 
elastically, but that a significant amount of fermion number can stick on the boundary. In 
this case, the K and B modes are anticorrelated in g, but in contrast to the kink-antikink 
collisions, very little fermion number is transfered to the KE mode or radiation. For the 
values of g employed here, low v favours transfer to the boundary, whereas for high v the 
fermions stay on the kink and are carried away from the boundary again. 

We in fact also solved for the evolution of the initial B, A, KE, and AE modes, and 
unsurprisingly found that fermions initially in these modes are also distributed on all modes 
after collision. In particular, a fermion localised on the boundary can be carried away by 
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a kink bouncing on this boundary. For brevity and to focus on our main aim, we did not 
carry out a detailed exposition of all mode combinations. It is however straightforward to 
do so. 

We found that the numerical implementation of the zizBC boundary conditions require 
some care, and that discretisation errors can be significant. We dealt with this by using 
iterative methods, higher order derivatives and rather fine lattices, Ddr = 0.0025. 

One could consider including more and more fermion modes in the spectrum, in order 
to track the "lost" fermions. In the end, this would lead to a full quantum (Hartree) 
treatment [^], and one would be able to include the back-reaction on the scalar self- 
consistently. For the purpose of this paper, we found the K, A, KE, AK and B modes to 
be sufficient. 

We started by arguing that our 1+1 model could be lifted to give results for a 4+1 
spacetime. While this is true, one has ignored any effects that depend on the directions 
within the brane. In particular we have not presented any results about the k dependence of 
the Bogoliubov coefficients (where k refers to the wave-number in the brane world- volume.) 
We hope to present results on this in a future publication. 

Another natural extension of this work is to consider more than one species of fermions 
and/or more scalars to include effects like C and CP violation in the scalar- fermion inter- 
action. This would for instance be relevant to electroweak baryogenesis, where at a first 
order phase transition, a domain wall sweeps through a plasma, reflecting fermions off it 
in a CP-violating way. 
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